Instability and dynamics of two nonlinearly coupled laser beams in a plasma 
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We investigate the nonlinear interaction between two laser beams in a plasma in the 
weakly nonlinear and relativistic regime. The evolution of the laser beams is governed by 
two nonlinear Schrodinger equations that are coupled with the slow plasma density response. 
We study the growth rates of the Raman forward and backward scattering instabilities as 
well of the Brillouin and sclf-focusing/modulational instabilities. The nonlinear evolution of 
the instabilities is investigated by means of direct simulations of the time-dependent system 
of nonlinear equations. 

PACS numbers: 52.35.Hr, 52.35.Mw, 52.38. Bv, 52.38. Hb 

I. INTRODUCTION 

The interaction between intense laser beams and plasmas leads to a variety of different insta- 



s to a variety or 

mm 



scattering and 



bilities, including Brillouin and Raman forward and backward 
modulational instabilities. In multiple dimensions we also have filamentation and side-scattering 
instabilities. Relativistic effects can then play an important role Q, ||3, Q]. When two laser beams 
interact in the plasma, we have a new set of phenomena. An interesting application is the beat- 
wave accelerator, in which two crossing beams with somewhat different frequencies can accelerate 
electrons to ultra-relativistic speeds via the ponderomotive force acting on the electrons. The 
modulational and filamentation instabilities of multiple co-propagating electromagnetic waves can 
be described by a system of coupled nonlinear Schrodinger equations from which the nonlinear 
wave coupling and the interaction between localized light wave packets can be easily studied 
Two co-propagating narrow laser beams may attract each other and spiral around each other [lCJ 
or merge Counter-propagating laser beams detuned by twice the plasma frequency can, at 



relativistic intensities, give rise to fast plasma waves via higher-order nonlinearities 



la lid 



li. At 



relativistic amplitudes, plasma waves can also be excited via beat wave excitation at frequencies 



2 



different from the electron plasma frequency, with applications to efficient wake-field accelerators 
[~L5|. The relativistic wakefield behind intense laser pulses is periodic in one-dimension .1.6] and 
shows a quasi-periodic behavior in multi-dimensional simulations [171 ] . Particle-in-cell simulations 
' lave demonstrated the generation of large-amplitude plasma wake-fields by colliding laser pulses 
lsl ] or by two co-propagating pulses where a long trailing pulse is modulated efficiently by the 
periodic plasma wake behind the first short pulse 

In the present paper, we consider the nonlinear interaction between two weakly relativistic 
crossing laser beams in plasmas. We derive a set of nonlinear mode coupled equations and nonlinear 
dispersion relations, which we analyze for Raman backward and forward scattering instabilities as 
well as for Brillouin and modulation/self- focusing instabilities. 

II. NONLINEAR MODEL EQUATIONS 

We consider the propagation of intense laser light in an electron-ion plasma. The slowly varying 
electron density perturbation is denoted by n es \. Thus, our starting point is the Maxwell equation 

« ™ ^< , ^ IdE 

V x B = — -(n + n es i)ev + - — . (1) 

The laser field is given in the radiation gauge, B = Vxi and E = —(l/c)dA/dt. Since dp e /dt = 
—eE, we thus have p e = eA/c. Moreover, p e = m e ^v e , where m e is the electron rest mass and 
7 = (1 — fg/c 2 ) -1 / 2 is the relativistic gamma factor, so that 

_ el / 2e 2 |l| 2 V V2 

! + -ir • ( 2 ) 
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For weakly relativistic particles, i.e. e \A\ /m e c <C 1, we can approximate @ by 

_ el ( e 2 \A\ 2 \ 
m e c \ m^c 1 J 

With these prerequisites, Eq. Q becomes 

" ^V 2 ) A + ul (1 + N s ) A-ufri^A = 0, (4) 

where u p o = (47rnoe 2 /m e ) 1//2 is the electron plasma frequency and we have denoted N s = n es i/riQ. 

Next, we divide the vector potential into two parts according to A = A\ + Ai , representing the 
two laser pulses. We also consider the case A\ ■ A<i ~ 0. With this, we obtain from Q) the two 
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coupled equations 



o2 \ 2 

° ,2v72\ A . , , .2 , at \ A ,.2 e (\ A |2 , I 7 _ |2\ 



_ _ c ^v 2 J A x + o^o (1 + 7V S ) A\ - ^-^(l^l 2 + I^IOAi = 0, (5a) 



and 

)2 



JL _ C 2 V 2^) A 2 + ^ (1 + ^ ^ _ ^ JL^p + l^p)^ = 0. (5b) 

Assuming that Aj is proportional to exp(ikj ■ r — iuijt) where ujj ^> \d/dt\, we obtain in the slowly 

varying envelope approximation two coupled nonlinear Schrodinger equations 

/ q \ _ _ _ e 2 --^ 

-2iwi ( — + • Vj Ai - c 2 V 2 Ai + w 2 iV s Ai - ^-^flAil 2 + |A 2 | 2 )Ai = 0, (6a) 

and 

-2w 2 ( — + v g2 • Vj A 2 - c 2 V 2 A 2 + co 2 p0 N s A 2 - ^-^flAxI 2 + |A 2 | 2 )A 2 = 0, (6b) 

where u 9 j = kjC 2 /ujj is the group velocity and = (a; 2 + (?k 2 ) 1 / 2 is the electromagnetic wave 
frequency. 

In order to close (JBJ), we next consider the slow plasma response. Here we may follow two 
routes. First, if we assume immobile ions, the slowly varying electron number density and velocity 
perturbations satisfy the equations 

dn es i 



at 

and 



+ n V • v es i = 0, (7) 



%i + -^VdAil 2 + |A 2 | 2 ) = -^V0 S - -^Vn esl (8) 
where T e is the electron temperature, together with the Poisson equation 

V 2 4> s = Series!. (9) 
Thus, combining ( |7| )— (|9" |) together with the vector potential decomposition, we obtain 

1 ° " .2 vt2 



\ dt 2 -re« ' m 2 r > 



34 e V 2 +o; 2 N s = V 2 (\A^ + |A 2 | 2 ), (10) 



where the electron thermal velocity is denoted by vxe = (Te/rrie) 1 ^ 2 . 

Second, if the electrons are treated as inertialess, we have in the quasi-neutral limit rtj s i = 
n esl = n s i 

V(|Ai| 2 + |A 2 | 2 ) =n eV<t> s -T e Vn sl , (11) 



m e c 2 



and 



dvi s i - 

n mi—^- = -n o eV0 s - 3TjVn s i. (12) 



Adding Eqs. (fTT|) and (fT2j). we obtain 

f)i) n p 2 - 

nom^ + ^V(|li| 2 + L4 2 | 2 ) + (T e + 3T,)Vn sl = 0, (13) 
ot m e c l 

which should be combined with 



at 

to obtain 



<) "- 1 +n o V-v isl = 0, (14) 



( S - ^) iV s = -^-,V 2 (|i 1 | 2 + L4 2 | 2 ), (15) 



\dt 2 s J S m e mi6 



where the sound speed is c s = y (T e + 3Tj)/mj and Tj is the ion temperature. 

III. COUPLED LASER BEAM AMPLITUDE MODULATION THEORY 

We shall consider, successively, Eqs. (6a, b) combined with ()1U|) (Case I: Raman scattering) or 
with p5[) (Case II: Brillouin scattering). 

A. Evolution equations 

Setting V — > i-K" and d/dt — > —iQ into the equations for the plasma density responses, we 
obtain 

iV s = ao(|i*i| 2 + |i' 2 | 2 ), (16) 

where, for Case I: 

e 2 K 2 

°0 = 2~~2 7v? o ts9 2 2~~ ' (17a) 

mic z U z - 2>K i vt^ e — a;p 

and for Case II: 
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The expressions (|16|) and (|17|) derived above provide the slow plasma response for any given 
pair of fields {Aj} (j = 1, 2). The latter now obey a set of coupled equations, which are obtained 
by substituting (|T6|) into ©, 

2iui (2. + v gl + c 2 V 2 A 1 + lo 2 p0 (-^ - a^j (|Ai| 2 + \A 2 \ 2 )A 1 = , (18a) 



and 



2*w 2 ^ + t; 92 -V^I 2 + c 2 V 2 A 2 +^ ^_|__ ao ^ (l^l 2 + |1 2 | 2 )1 2 = 0, (18b) 



For convenience, Eqs. (|18b) and (|18b) are cast into the reduced form as 

2»wi + u 5 i • V ] A x + c 2 V 2 A 1 + Q (|li| 2 + |A 2 | 2 )Ai = , (19a) 



and 



2iuo 2 ( ^ + % 2 • Vj A 2 + c 2 V 2 I 2 + Q (|li| 2 + |I 2 | 2 )A 2 = , (19b) 

where Aj has been normalized by m e c 2 /e and where the nonlinearity /coupling coefficients are 

K 2 c 2 

' J Te ~ ^ P 



and 



Q = ^o U-^T^ro , (20b) 



m e K 2 c 2 

J p°\'~ n 2 - K 2 c 2 J ' 
for stimulated Raman (Case I) and Brillouin (Case II) scattering, respectively. We observe that 
the expressions (|2Ub ) and (J^UJd) may be either positive or negative, depending on the frequency f2, 
prescribing either the modulational instability or the Raman and Brillouin scattering instabilities 

The two nonlinear wave equations are identical upon an index (1,2) interchange, and coincide 
for equal frequencies u>i = uj 2 . 



B. Nonlinear dispersion relation 

We now investigate the parametric instabilities of the system of equations (|19a|) and ()19b|) . 
Fourier decomposing the system by the ansatz Aj = [Ajq + Aj + exp(iK ■ r — i£lt) + Aj- exp(—iK ■ 
r + iSlt))} exp(— itlot), where |A/o| 3> anc ^ sorting for different powers of exp(ii^ • r — i£lt), 

we find the nonlinear frequency shift 

%o = -Qx=o(|^io| 2 + \A 20 \ 2 )/2u; v (21) 
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where Qk=o denotes the expression for Q with K = 0. For the nonlinear wave couplings, we have 
from the system of equations 

D 1+ X 1+ + Q\A 10 \ 2 (X 1+ + X i+ + X 2+ + X 2 _) = 0, (22a) 

L>i_Xi_ + Q|^i | 2 (X 1+ + Xi+ + X 2+ + X 2 -) = 0, (22b) 

D 2+ X 2+ + Q\A 20 \ 2 (X 1+ + X 1+ + X 2+ + X 2 _) = 0, (22c) 

Di-Xt- + Q|^2o| 2 (^i+ + Xi+ + X 2+ + X 2 _) = 0, (22d) 

where the unknowns are X\ + = A\ -A\ + , X\- = A\q-A\_, X 2+ = A 20 -A 2+ , and X 2 - = A 2 q-A 2 _. 
The sidebands are characterized by 

Dj± = ±2[ujjn - c 2 kj ■ K] - c 2 K 2 , (23) 

where we have used v g j = ckj/ujj. The solution of the system of equations (|22|) yields the nonlinear 
dispersion relation 

which relates the complex-valued frequency O to the wavenumber K. Equation (|24|) covers Raman 
forward and backscattering instabilities, as well as the Brillouin backscattering instability or the 
modulational/self- focusing instability, depending on the two expressions for the coupling constant 
Q. If either \A\o\ or IA20I is zero, then we recover the usual expressions for a single laser beam in 
a laboratory plasma, or for a high-frequency radio beam in the ionosphere 21]. 



IV. NUMERICAL RESULTS 



We have solved the nonlinear dispersion relation (|24l) and presented the numerical results in Figs. 
1-5. In all cases, we have used the normalized weakly relativistic pump wave amplitudes Ajq = 0.1 
with different sets of wavenumbers for the two beams. The nonlinear couplings between the laser 
beams and the Langmuir waves, giving rise to the Raman scattering instabilities (Case I), are 
considered in Figs. 1 and 2. The instability essentially obeys the matching conditions ujj = lo s + 0, 
and kj = k s + K, where ujj and kj are the frequency and wavenumbers of the pump wave, uj s and k s 
are the frequency and wavenumbers for the scattered and frequency downshifted electromagnetic 
daughter wave, f2 and K are the frequencies of the Langmuir waves, and where the light waves 
approximately obey the linear dispersion relation, ujj = (uj 2 + fcjc 2 ) 1 / 2 , uj s = (uj 2 + k^c 2 ) 1 / 2 and 
the low-frequency waves obey the Langmuir dispersion relation = (uj 2 + 3-fT 2 Wy e ) 1//2 . We thus 
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have the matching condition (u; 2 + kjc 2 ) 1 / 2 = [w 2 + (kj - K) 2 c 2 ]^ 2 + (u; 2 + SK 2 ^^) 1 / 2 , which 
in two-dimensions relates the components K y and K z of the Langmuir waves to each other, and 
which gives rise to almost circular regions of instability, as seen in Figs. 1 and 2. In the upper left 
and right panels of Fig. 1, we have assumed that the single beams A\ and Ai propagate in the 
y and z direction, respectively, having the wavenumber (ki y ,ki z ) = (6, 0) and (A^y, k2 Z ) = (0, 4), 
respectively. We can clearly see a backward Raman instability, which for the beams A\ and A2 
have maximum growth rates at (K y ,K z ) = (2k\ y , 0) = (12, 0)uj p q/c and (K y , K z ) = (0, 2k2 Z ) = 
(0, 8)uj p o/c, respectively. The backward Raman instability is connected via the obliquely growing 
wave modes to the forward Raman scattering instability that has a maximum growth rate (much 
smaller than that of the backward Raman scattering instability) at the wave number K ss uj pe /c in 
the same directions as the laser beams. In the lower panels, we consider the two beams propagating 
simultaneously in the plasma, at a right angle to each other (lower left panel) and in opposite 
directions (lower right panel) . We see that the dispersion relation predicts a rather weak interaction 
between the two laser beams, where the lower left panel shows more or less a superposition of the 
growth rates in the two upper panels. The case of two counter-propagating laser beams (lower right 
panel) also shows a weak interaction between the two beams. For the case of equal wavelengths of 
the two pump waves, as shown in Fig. 2, we have a similar scenario as in Fig. 1. The lower left 
panel of Fig. 2 shows that the growth rate of two interacting laser beams propagating at a right 
angle to each other is almost a superposition of the growth rates of the single laser beams displayed 
in the upper panels of Fig. 2. Only for the counter-propagating laser beams in the lower right 
panel we see that the instability regions have split into broader and narrower bands of instability, 
while the magnitude of the instability is the same as for the single beam cases. 

We next turn to the Brillouin scattering scenario (Case II), in which the laser wave is scattered 
against ion acoustic waves, displayed in Figs. 3 and 4. In the weakly nonlinear case, we have three- 
wave couplings in the same manner as for the interaction with Langmuir waves, and we see in both 
Figs. 3 and 4 that the instability has a maximum growth rate in a narrow, almost circular band 
in the (K y , K z ) plane. In the upper two panels, we also see the backscattered Brillouin instability 
with a maximum growth rate at approximately twice the pump wavenumbers, but we do not have 
the forward scattered instability. Instead, we see a broadband weak instability in all directions and 
also perpendicular to the pump wavenumbers. A careful study shows that the perpendicular waves 
are purely growing, i.e. there may be density channels created along the propagation direction of 
the laser beam. In the lower panels of Figs. 3 and 4, we display the cases with interacting laser 
beams. Also in the case of Brillouin scattering, the nonlinear dispersion relation predicts a rather 
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weak interaction between the two beams, where the instability regions of the two beams are more 
or less superimposed without dramatic differences in the growth rates. 

In order to investigate the nonlinear dynamics of the interacting laser beams in plasmas, we have 
carried out numerical simulations of the reduced system of equations © in two spatial dimensions, 
and have presented the results in Figs. 5-8. In these simulations, we have used as an initial 
condition that either A\ has a constant amplitude of 0.1 and A2 has a zero amplitude, or that both 
beams have a constant amplitude of 0.1 and that they initially have group velocities at a right angle 
to each other. Due to symmetry reasons, it is sufficient to simulate one vector component of Aj, 
which we will denote Aj (j = 1,2). The background plasma density is slightly perturbed with a low- 
level noise (random numbers). We first consider stimulated Raman scattering, displayed in Figs. 
5 and 6. The single beam case in Fig. 5 shows a growth of density waves mainly in the direction of 
the beam, while a standing wave pattern is created in the amplitude of the electromagnetic wave 
envelope, where maxima in the laser beam amplitude is (roughly) correlated with minima in the 
electron density. This is in line with the standard Raman backscattering instability. The simulation 
is ended when the plasma density fluctuations are large and self-nonlinearities and kinetic effects 
are likely to become important. In Fig. 6, we show the case with the two beams crossing each 
other at a right angle. In this case, the wave pattern becomes slightly more complicated with 
local maxima of the laser beam envelope amplitude correlated with local minima of the electron 
density. However, this pattern is very regular and there is no clear sign of nonlinear structures 
in the numerical solution. We next turn to the case of stimulated Brillouin scattering, presented 
in Figs. 7 and 8. In this case, the waves grow not only in the direction of the laser beam but 
also, with almost the same growth rate, obliquely to the propagation direction of the laser beam. 
We see in the single beam case, presented in Fig. 7, that the envelope of the ion beam becomes 
modulated in localized areas both in y and z directions, and in the nonlinear phase at the end 
of the simulation, the laser beam envelope has local maxima correlated with local minima of the 
ion density. For the case of two crossed laser beams, displayed in Fig. 8, we see a more irregular 
structure of the instability and that at the final stage, local "hot spots" are created in which large 
amplitude laser beam envelopes are correlated with local depletions of the ion density. 



V. SUMMARY 



In summary, we have investigated the instability and dynamics of two nonlinearly interact- 
ing intense laser beams in an unmagnetized plasma. Our analytical and numerical results reveal 
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that stimulated Raman forward and backward scattering instabilities are the dominating nonlinear 
processes that determine the stability of intense laser beams in plasmas, where relativistic mass 
increases and the radiation pressure effects play a dominant role. Our nonlinear dispersion relation 
for two interacting laser beams with different wavenumbers predicts a superposition of the instabil- 
ities for the single beams. The numerical simulation of the coupled nonlinear Schrodinger equations 
for the laser beams and the governing equations for the slow plasma density perturbations in the 
presence of the radiation pressures, reveal that in the case of stimulated Raman scattering, the 
nonlinear interaction between the two beams is weaker than for the case of stimulated Brillouin 
scattering. The latter case lead to local density cavities correlated with maxima in the electro- 
magnetic wave envelope. The present results should be useful for understanding the nonlinear 
propagation of two nonlinear ly interacting laser beams in plasmas, as well as for the acceleration of 
electrons by high gradient electrostatic fields that are created due to stimulated Raman scattering 
instabilities in laser-plasma interactions. 



10 



P. K. Shukla, N. N. Rao, M. Y. Yu, and N. L. Tsintsadze, Phys. Rep. 135, 1 (1986). 
A. Sjolund and L. Stenflo, Appl. Phys. Lett. 10, 201 (1967). 

M. Y. Yu, K. H. Spatschck, and P. K. Shukla, Z. Naturforsch. A 29, 1736 (1974). 

P. K. Shukla, M. Y. Yu, and K. H. Spatschck, Phys. Fluids 18, 265 (1975). 

P. K. Shukla and L. Stenflo, Phys. Rev. A 30, 2110 (1984). 

N. L. Tsintsadze and L. Stenflo, Phys. Lett. A 48, 399 (1974). 

C. E. Max, J. Arons, and A. B. Langdon, Phys. Rev. Lett. 33, 209 (1974). 

P. K. Shukla, Phys. Scripta 45, 618 (1992). 

L. Berge, Phys. Rev. E 58, 6606 (1998). 

C. Ren, B. J. Duda, and W. B. Mori, Phys. Rev. E 64, 067401 (2001). 
Q.-L. Dong, Z.-M. Sheng, and J. Zhang, Phys. Rev. E 66, 027402 (2002). 
M. N. Rosenbluth and C. S. Liu, Phys. Rev. Lett. 29, 701 (1972). 
G. Shvcts and N. J. Fisch, Phys. Rev. Lett. 86, 3328 (2001). 

R. Bingham, J. T. Mendonca and P. K. Shukla, Plasma Phys. Control. Fusion 46, Rl (2004). 
G. Shvcts, Phys. Rev. Lett. 93, 195001 (2004). 

V. I. Bcrczhiani and I. G. Murusidzc, Phys. Lett. A 148, 338 (1990). 

F. S. Tsung, R. Narang, W. B. Mori, R. A. Fonseca, and L. O. Silva, Phys. Rev. Lett. 93, 185002 
(2004). 

K. Nagashima, J. Koga, and M. Kando, Phys. Rev. Lett. 64, 066403 (2001). 

Z.-M. Sheng, K. Mima, Y. Setoku, K. Nishihara, and J. Zhang, Phys. Plasmas 9, 3147 (2002). 

N. L. Tsintsadze, D. D. Tskhakaya, and L. Stenflo, Phys. Lett. A 72, 115 (1979). 

L. Stenflo, Phys. Scripta T30, 166 (1990); ibid T107, 262 (2004). 



11 



Figure captions 

FIG. 1: The normalized (by uj p q) growth rates due to stimulated Raman scattering (Case I) 
for single laser beams (upper panels) and for two laser beams (lower panel), as a function of the 
wave numbers K y and K z . The upper left and right panels show the growth rate for beam A\ 
and A-2,., respectively, where the wave vector for A\ is (k y , k z ) = (6, 0)uj p o/c and the one for A2 
is (k y , k z ) = (0, 4)w p o/c, i-e. the two beams are launched in the y and z directions, respectively. 
In the lower left panel, A\ and A2 are launched simultaneously at a perpendicular angle to each 
other, and in the lower right panel, the two beams are counter-propagating. We used the normalized 
amplitudes |Alo| = | A20 1 = 0.1 and the electron thermal speed vr e = 0.01c. 

FIG. 2: The normalized (by to p o) growth rates due to stimulated Raman scattering (Case I) 
for single laser beams (upper panels) and for two laser beams (lower panel), as a function of the 
wave numbers K y and K z . The upper left and right panels show the growth rate for beam A\ 
and A2, respectively, where the wavenumber for A\ is (k y , k z ) = (5, 0) uj p q/c and the one for A2 is 
(k y , k z ) = (0, 5)lo p q/c. In the lower left panel, two beams are launched at a perpendicular angle 
to each other, and in the lower right panel, the two beams are counter-propagating. We used the 
normalized amplitudes \Aiq\ = |^2o| = 0.1 and the electron thermal speed vx e = 0.01c. 

FIG. 3: The normalized (by oj p q) growth rates due to stimulated Brillouin scattering (Case II) 
for single laser beams (upper panels) and for two laser beams (lower panel), as a function of the 
wave numbers K y and K z . The upper left and right panels show the growth rate for the beam A\ 
and A2, respectively, where the wave number for A\ is (k y , k z ) = (6, Q)lv p q/c and the one for A2 
is (k y , k z ) = (0, 4)uj p o/c. In the lower left panel, two beams are launched at a perpendicular angle 
to each other, and in the lower right panel, the two beams are counter-propagating. We used the 
normalized amplitudes | A10 1 = I-A20I = 0.1, the ion to electron mass ratio mi/m e = 73440 (Argon), 
and the ion sound speed c s = 3.4 x 10 -5 c. 

FIG. 4: The normalized (by uj p q) growth rates due to stimulated Brillouin scattering (Case II) 
for single laser beams (upper panels) and for two laser beams (lower panel), as a function of the 
wave numbers K y and K z . The upper left and right panels show the growth rate for beam A\ 
and A2, respectively, where the wavenumber for A\ is (k y , k z ) = (5, 0) oj p q/c and the one for A2 is 
(k y , k z ) = (0, 5)uj p o/c. In the lower left panel, two beams are launched at a perpendicular angle 
to each other, and in the lower right panel, the two beams are counter-propagating. We used the 
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normalized amplitudes | Aio | = I-A20I = 0.1, the ion to electron mass ratio mi/m e = 73440 (Argon), 
and the ion sound speed c s = 3.4 x 10~ 5 c. 

FIG. 5: The amplitude of a single laser beam (left panels) and the electron density N s 
(right panels) involving stimulated Raman scattering (Case I), at times t = 1.0 uj~q , t = 30 ^"q 1 
and t = OOWpQ 1 (upper to lower panels). The laser beam initially has the amplitude A\ = 0.1 and 
wavenumber {k\ y , k\ z ) = (0, 5)uj p q/c. The electron density is initially perturbed with a small- 
amplitude noise (random numbers) of order 10~ 4 . 

FIG. 6: The amplitude of two crossed laser beams, \A\ = (\Ai\ 2 + | ^4 2 1 2 ) (left panels) and 
the electron density N s (right panels) involving stimulated Raman scattering (Case I), at times 
t = l.OcjpQ 1 , t = 30ujpQ and t = 60 (upper to lower panels). The laser beams initially have the 
amplitude A\ = A2 = 0.1, and A\ initially has the wavenumber (k± y , k± z ) = (0, 5)uj p q/c while A2 
has the wavenumber (A^y, ^2z) = (5, 0)w p o/c The electron density is initially perturbed with a 
small-amplitude noise (random numbers) of order 10 -4 . 

FIG. 7: The amplitude of a single laser beam \A±\ (left panels) and the electron density N s 
(right panels) involving stimulated Brillouin scattering (Case II), at times t = 1.5lo~q, t = 600 uj~q 
and t = 1200 WpQ 1 (upper to lower panels). The laser beam initially has the amplitude A\ =0.1 and 
wavenumber {k\ y , k± z ) = (0, 5) uj p q/c. The ion density is initially perturbed with a small-amplitude 
noise (random numbers) of order 10 -4 . 

FIG. 8: The amplitude of two crossed laser beams, \A\ = (\Ai\ 2 + | ^4 2 1 2 ) 1 / 2 (left panels) and 
the electron density N s (right panels) involving stimulated Brillouin scattering (Case II), at times 
t = l.OcjpQ 1 , t = SOujpQ and t = 60 a;^ 1 (upper to lower panels). The laser beams initially have the 
amplitude A\ = A2 = 0.1, and A\ initially has the wavenumber (ki y , k\ z ) = (0, 5)iv p o/c while A2 
has the wavenumber (k2 y , k2 Z ) = (5, 0)uj p q/c. The electron density is initially perturbed with a 
small-amplitude noise (random numbers) of order 10~ 4 . 
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